On the equivalence of pairing correlations and intrinsic vortical currents 

in rotating nuclei 
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The present paper establishes a hnk between pairing correlations in rotating nuclei and collective 
vortical modes in the intrinsic frame. We show that the latter can be embodied by a simple S- 
type coupling a la Chandrasekhar between rotational and intrinsic vortical collective modes. This 
results from a comparison between the solutions of microscopic calculations within the HFB and 
the EIF Routhian formalisms. The HF Routhian solutions are constrained to have the same Kelvin 
circulation expectation value as the HFB ones. It is shown in several mass regions, pairing regimes, 
and for various spin values that this procedure yields moments of inertia, angular velocities, and 
current distributions which are very similar within both formalisms. We finally present perspectives 
for further studies. 
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I. INTRODUCTION AND METHODS 

In a recent paper [1] , some of the authors of the present 
article conjectured that the well-known dynamical effects 
of pairing correlations in rotating nuclei, explicited in the 
significant decrease of the moments of inertia from rigid 
body values [2,3], could be pictured as a vortical intrinsic 
collective motion coupled to the global rotation of the so- 
called intrinsic (body-fixed) frame relatively to the lab- 
oratory system. Namely, one intuitively would think of 
such a mode as aligned and counter-rotating with respect 
to the global angular velocity so as to produce the above 
quoted reduction of the moments of inertia. Further- 
more, if one wants to preserve the nuclear shape in the 
presence of this intrinsic vortical mode the latter should 
be tangential to the grossly-defined nuclear surface. 

As often suggested after R. Y. Cusson [4] (see also Ref. 
[5] for some relevant references) , a simple ansatz for such 
a coupling is provided by the S-ellipsoid fluid dynamics 
as studied classically in great details by Chandrasekhar 
in the context of celestial self-gravitating objects [6]. It 
makes use of a velocity field which is linear in the coordi- 
nates and produces a coupling of a global rotation with 
an intrinsic mode corresponding to a motion tangential 
to the body surface which is supposed to be eUipsoidal. 
Furthermore, the vorticities (i.e., the curl) of the two ve- 
locity fields are assumed to be aligned or anti-aligned 
along a principal axis of the ellipsoid. As shown, e.g., in 
Ref. [7] (see Fig. 1 there) such a simple parameterization 
is capable of yielding a variety of modes in the laboratory 
frame from rigid-body rotation to irrotational modes as 
well as various shear modes. 



The connection of this S-ellipsoid fluid dynamics with 
pairing correlations may be discussed at various levels. 
First, one may note that a tangential intrinsic mode 
which is not changing the body shape (in the r-sector 
of the phase-space) is only able to redistribute the den- 
sity function in the p-sector of the phase-space. Clearly, 
up to self-consistent couplings between the momentum 
redistribution and the body shape, the BCS pairing cor- 
relations behave exactly in the same way (redistributing 
the single-particle momenta). Second, in the semiclassi- 
cal (Thomas- Fermi) calculations of Ref. [8] , upon widely 
varying the pairing strength in some deformed rare-earth 
rotating nuclei, the authors have unambiguously exhib- 
ited current patterns very similar indeed to the classical 
current patterns of S-ellipsoids. Such currents had been 
shown before to exist in paired solutions at finite spin 
[9,10]. They were however embedded in the midst of 
shell-effect — generated intrinsic currents whose existence 
had been found long ago (see the non-paired rotating so- 
lutions of [11]) and which will be further discussed in 
Section II A below. 

However, the arguments above sketched are purely 
qualitative. It is the aim of this paper to investigate 
how quantitatively correct they might be. 

To achieve this goal, it is first necessary to find a tool 
to quantify the intrinsic vortical content of a given cur- 
rent distribution. Within the context of S-ellipsoid fluids, 
it is rather natural to use the Kelvin circulation whose 
components are defined by (see, e.g., [12]) 
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where e^fc is the totally antisymmetrical third rank ten- 
sor and the aj 's are characteristic lengths associated with 
the ellipsoidal shape of the nuclei (e.g., the semi-axes). 
These lengths arc proportional to \/{xf)- These op- 
erators clearly appear as the components of a doubly- 
strctchcd (in both r and p) orbital angular momentum. 
Physically, they correspond to a rotation in the intrinsic 
frame of a sphere obtained by stretching the ellipsoidal 
intrinsic distribution. Indeed as shown by Lebowitz (see, 
e.g., Ref. [6]), the S-ellipsoid coupling mode may be un- 
derstood as a global rotation along some principal axis 
71, followed by the stretching to a sphere and the rotation 
(according to the same axis u of the considered ellipsoid) 
as above described, and finally by the inverse stretching 
back to the original ellipsoidal shape. 

To check the validity of the above mentioned con- 
jecture, we compare the results of two types of micro- 
scopic calculations, one involving pairing correlations in 
a purely rotating formalism, the other with no pairing 
but imposing the S-ellipsoid dynamics for the global ro- 
tation and intrinsic vortical mode velocity fields as well 
as for their coupling. 

The pairing correlations htc treated a, la Hartrec-Fock- 
Bogoliubov (HFB) using a Skyrme effective interaction 
in the particle-hole channel and a seniority force in the 
pairing channel within a code developed by H. Laftchicv 
[13,14]. The triaxial character of the solutions arising 
from the time-reversal symmetry breaking inherent to 
rotating nuclei is treated within the Fourier decomposi- 
tion formalism described in Ref. [15]. Namely, the wave- 
functions and relevant densities are defined in cylindri- 
cal coordinates on a mesh in the two spatial directions 
(r, z) and as Fourier series in the azimuthal angle 9. 
An approximate description of particle-number projec- 
tion within the Lipkin-Nogami (LN) scheme [16] is also 
available in the code. 

The S-cUipsoid coupling case is described within the 
Hartree-Fock-Skyrme generalized Routhian formalism of 
Refs. [15,17]. This formalism makes it possible to de- 
scribe the S-cUipsoid dynamics upon using a double con- 
straint on both angular momentum and Kelvin circula- 
tion via two angular velocities 17 and w (relevant to each 
collective mode) along the first-axis. 

These two microscopic approaches have been used in 
the following scheme. In a first step, we have used the 
HFB formalism with a constraint on angular momentum 
using realistic pairing strengths. We thus generate vorti- 
cal currents in the intrinsic frame which can be quantified 
for instance by the mean value of the Kelvin circulation 
operator measured in the intrinsic frame and defined as 



{K,) = I (f xj)i d^r. 



(2) 



There r is the stretched position vector and j is the 

stretched current in the intrinsic frame whose compo- 
nents are obtained after subtracting the rotating-part to 
the usual current-density j as 



TABLE I. Pairing strength G„ and Gp for neutrons and 
protons as well as energy cutoff n (in MeV) for each of the 

three studied nuclei. The matrix elements gq of the seniority 
force are defined by the usual prescription Qq — Gq/{11 + Nq). 
The configuration space includes all single-particle states 
whose energies lies below e + n where e is the chemical poten- 
tial. 
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^^^Hg 


2B4No 


Gn 


14.4 


11.5 


13.8 


Gp 


14.4 


11.5 


15.0 


Cutoff (yU) 


5 


10 


7 
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where j is the usual quantum mechanical current density 
defined as 
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The mean value / of the angular momentum operator is 
defined through 



(5) 



where p is the spin- vector density. 

In a second step, we perform pure Hartree-Fock cal- 
culations (i.e., without pairing) with a double constraint 
on both angular momentum and Kelvin circulation to 
obtain the same values / and (ATi) as in the preceding 
step. In the following, this type of calculations will be 
referred to as HF+V calculations. 

The paper will be organized as follows. In the next Sec- 
tion we will present and compare the results of the numer- 
ical calculations yielded by the two formalisms, namely 
HFB and HF-I-V in several mass regions. Section III will 
be (kn'oted to a summary of the salient features exhibited 
in Section II and to some perspectives for further studies. 



II. RESULTS 

Numerical calculations have been performed for three 
heavy nuclei: two in the A ~ 150 and A ~ 190 superde- 
formation regions (namely for the ^^'^Gd and the ^^^Hg 
nuclei) including LN approximate particle-number pro- 
jection scheme and one very heavy nucleus (^^'^No) with- 
out LN projection. The SkM* parameterization [18] of 
the Skyrme effective force has been used. The pairing 
strengths Gq for the charge state q arc given for each cal- 
culated nuclei in table I together with the corresponding 
single-particle configuration space in use. For all three 
nuclei, kinetic (3^^^) and dynamical (3^^^) moments of 
inertia are calculated using the formulae 
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FIG. 1. Dynamical (upper panel) and kinetic (lower panel) moments of inertia (in ft^MeV'^) for the three considered nuclei 
as functions of the angular velocity H in MeV. The conventions in use are the following: HF value (dotted line), HFB value 
(full line), and HF+V value (dashed line). Experimental data are represented as full circles, excepted for the ^^"Gd kinetic 
moment of inertia with the assumption /init = 34 as opened circles (lower left panel). The Lipkin-Nogami correction has been 
applied for ^^OGd and ^''^Hg. 
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As was already mentioned in a previous paper [14], us- 
ing these expressions for the moments of inertia, we as- 
sume that they can be computed as well from the energy- 
corrected LN wavcfunctions. Wc also assume that the 
angular velocity dependence of the estimate of LN A2 
parameter is small, otherwise 3^^^ would have to be eval- 
uated through second derivatives of the energy with re- 
spect to fi. 



A. The yrast superdeformed band of ^ Gd 

Up to now, fourteen superdeformed bands have been 
experimentally found in ^^°Gd [19,20]. We will deal 
here with the yrast superdeformed band only. Using the 
Cranked-HF approach, it was impossible [15] to repro- 
duce the trends of the experimental kinetic and dynam- 
ical moments of inertia of this band as can be seen in 
Fig. 1 (dotted line). 

The results presented in the upper panel of Fig. 1 for 
the dynamical moment of inertia of this nucleus show a 
much better agreement with experimental data within 
the HFB+LN formalism as compared to the HF one. In- 
deed, we now reproduce the decrease of this moment of 



inertia along the band with theoretical values which are 
rather close to the experimental ones as was previously 
shown within pure HFB calculations by Bonche et al. 
[21]. In the lower panel, the agreement of HFB+LN ki- 
netic moment of inertia with the experimental data of 
Ref. [19], though much better than the HF one, is seen to 
be rather poor. However, the experimental value of this 
moment of inertia relies on an assumption of the lowest 
spin value of the observed band [see eq. (6)] which has 
been assumed in [19] to be /init = 32 h. Taking now for 
the initial spin of the band /init = 34 ?i shifts up the ex- 
perimental data and then yields a very good agreement 
with HFB+LN values of the moment of inertia. 

Let us now focus on the intrinsic vortical content of the 
current density. In the upper part of Fig. 2 the Kelvin cir- 
culation mean value is plotted as a function of the angular 
velocity for ^^°Gd. One may first notice that the Kelvin 
circulation mean value in the HF case is not zero, even 
though the intrinsic vorticity constraint is absent. It is 
negative (with a rather small absolute value of about 1 h) 
for the protons while it is positive and larger (about 2 h) 
for the neutrons. These contributions should most likely 
be attributed to small inhomogeneities in the density due 
to shell effects. It is almost constant for the neutrons and 
slightly decreasing with Q for the protons. Dubbing this 
current patterns as shell effects may also be understood 
in rather general terms. In Ref. [22] within purely semi- 
classical microscopic calculations, it has been shown that 
for finite systems small surface-peaked counter-rotating 
intrinsic currents do appear. This effect is analogous to 
the diamagnetism which is observed in an electronic sys- 
tem submitted to a constant magnetic field. Thence the 
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FIG. 2. Kelvin circulation mean value in ft units (upper panel) and pairing energy in MeV (lower panel) for the three 
considered nuclei as functions of the angular velocity f2 in MeV. The conventions in use are the following: proton contributions 
(dashed line), neutron contributions (dotted line), and total contributions (full line); in the upper panel: HF value (no symbols), 
HF+V value (opened circles), HFB value (full circles). The Lipkin-Nogami correction has been applied for ^^°Gd and ^^^Hg. 



origin of all intrinsic currents (in Routhian HF calcula- 
tions) which are not surface-peaked and counter-rotating 
(as obtained years ago in, e.g., [11]) are clearly and prop- 
erly to be attributed to shell effects. 

Comparing now the HF curves with the HFB ones, 
wc sec that the introduction of pairing correlations shifts 
down the Kelvin circulation mean value by about 2 h for 
each charge state at low angular velocities and still 0.5 H 
units at the end of the band. This decrease of the Kelvin 
circulation in the HFB case is related, as it has been 
already demonstrated in [1] , to the fact that the collective 
effect of pairing correlations is equivalent to a counter- 
rotating intrinsic vortical motion (see also Ref. [8]). The 
combination of the above mentioned HF variation of (Ki) 
with fl for protons and neutrons together with the regular 
decrease of the pairing correction as described previously 
yields an increase of the neutron contribution and an 
almost constant pattern for the protons. As a result the 
total Kelvin circulation exhibits an fl variation similar to 
the neutron one. 

This decreasing pairing effect on the ^^"Gd Kelvin cir- 
culation is nicely correlated, as can be seen on the lower 
panel of Fig. 2, with the decrease of the pairing energy 
(defined as a trace of the product of the abnormal den- 
sity with the gap potential). The latter is a well-known 
behavior, known as the Mottelson-Valatin effect [23], sim- 
ilar to the effect of a magnetic field on a superconductor 
below the critical point [24]. Indeed as the pairing en- 
ergy almost vanishes at high angular velocity, the Kelvin 
circulation in the HFB formalism almost reaches its HF 
value. 

As said in the Introduction, the Kelvin circulation is a 
powerful tool to measure the intrinsic vortical content of 



the current distribution. In Fig. 3 (^^°Gd bottom part) 
we present the current patterns in the intrinsic frame of 
two rotational states of the HFB supcrdcformed yrast 
band of Gd, namely for / = 32 and / = 46 h 
(corresponding to angular velocities O = 381 keV and 
r2 = 516 keV). In both cases, the neutron current pat- 
terns are rather disordered (due to shell structure effects). 
This is related to positive values of the neutron contribu- 
tion to the Kelvin circulations. On the contrary, proton 
current patterns are rather well oriented, at least around 
the nuclear surface, along elliptic lines parallel to the sur- 
face, the more so in the I = 32 H case where total Kelvin 
circulation and pairing energies are greater (in absolute 
value) than for I ~ 46 h. This trend is due to the fact 
(observed in Fig. 2) that the proton HF contribution to 
I {Ki) \ is very small. In the HF case (upper part), the cur- 
rent patterns show no particular order for both spins, and 
for protons as well as for neutrons. This clearly shows 
that the inclusion of pairing correlations favors the ap- 
pearance of S-ellipsoid-like currents in rotating nuclei. 

We then have performed HF+V calculations, namely 
constraining HF solutions to have / and (Ki) values [see 
eqs. (2) and (5)] identical to their HFB counterparts. 
The constraint in use is of an isoscalar type, that is con- 
straining on the total (neutron -I- proton) angular mo- 
mentum and Kelvin circulation. However, it can be seen 
in Fig. 2 that the proton and neutron contributions to the 
Kelvin circulation in HF+V calculations are very close 
to what is obtained in HFB calculations even though 
these contributions have not been constrained separately. 
Moreover, the current patterns for the HFB and HF+V 
approaches (see Fig. 3) are similar. Hence we have suc- 
cessfully grafted the superfluid part of the HFB currents 
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I=32h I = 46h I=lOh I = 2Qh I = 2h I = 2Qh 

FIG. 3. Current patterns in the intrinsic frame for two rotational states in the three considered nuclei. For each nucleus. 
The same arbitrary units for the length of \j\ are used for both spins. In all ellipses, the proton (resp. neutron) current patterns 
are represented in the top-right (resp. bottom-left) quarter. The global rotation direction is counter-clockwise. The formalisms 
used are HF, HF-I-V and HFB from top to bottom. The Lipkin-Nogami correction has been applied for ^^°Gd and ^^^Hg. 



into the HF calculations constraining merely the expec- 
tation values of the Ki isoscalar operator. 

The angular velocities within the three studied ap- 
proaches, namely HF, HFB and HF-I-V, are plotted in 
Fig. 4 against the angular momentum. It is seen that the 
HFB and HF-^V values fully overlap: 

^HFB — ^HF+V ■ (8) 

It is worth noting that these two angular velocities are 
greater than ^hf- This is not surprising, since the 
counter-rotating collective effect of pairing has to be com- 
pensated in HFB calculations at a given angular momen- 
tum by an increase of the angular velocity. However, it 
is striking to see that the identity of the two momenta of 
the current-distribution, namely / and (-ft'i) in HFB and 
HF-I-V calculations, leads to identical rotational angular 
velocities Q,. It proves indeed, more quantitatively than 
the plots of Fig. 3, that the current patterns in the two 
approaches are similar everywhere, which could hardly be 
demonstrated by merely looking at integrated quantities. 

The identity of the HFB and HF-I-V angular velocities 
seen in Fig. 4 yields an identity of the kinetic moments 
of inertia clearly demonstrated in eq. (6). As for the 
HF-hV dynamical moments of inertia, even though they 
do not fully coincide with their HFB counterparts, they 
are indeed very close and even reproduce better the ex- 
perimental data. In this nuclei, the HF-I-V and HFB 
formalisms yield solutions whose axial quadrupole mo- 
ments remain constant over the entire band (differing by 
less than 2 percent around Q2Q — 3900 fm-^). Hence the 
equivalence of the dynamical moments of inertia in the 



two formalisms is clearly free from any deformation effect 
and reveals similar superfluid properties. 



B. The yrast superdeformed band of Hg 

We have used also the HF, HFB+LN and HF+V for- 
malisms in order to describe the yrast superdeformed 
band of ^^^Hg and compare their results with the avail- 
able experimental data [25]. 

As already noted by Gall et al. [26] the very low 
amount of pairing correlations attained above / ~ 30 fi. 
(as exemplified by the pairing energies for both neutrons 
and protons in the corresponding lower part of Fig. 2) 
makes it necessary to correct for the particle number 
symmetry breaking. This has been performed upon using 
the approximate Lipkin-Nogami particle number projec- 
tion on top of the HFB calculation which has been found 
necessary to approach a proper description of the rota- 
tional band of this nuclei and avoid spurious drops in 
the moments of inertia due to sudden pairing-energy dis- 
appearance both for neutrons and protons. Whereas the 
LN prescription is known to be efficient in restoring some 
pairing correlations in weak-pairing regime as would be 
the case in HFB calculation for these nuclear states, the 
question of knowing whether or not it provides the correct 
amount of pairing correlations is still controversial (see, 
e.g., the discussion of Peru [27]). A more exact treatment 
of pairing correlations could be obtained for instance by 
generalizing to time-reversal symmetry breaking systems 
the method of Fillet et al. [28,29] which explicitly con- 



5 




/ (h) I (R) I (H) 

FIG. 4. Angular velocity fl (in MeV) as a function of the nuclear spin I (in h units) for the three considered nuclei. The 
conventions in use art; the following; HF value (dotted line), HFB value (full line), and HF+V value (dashed line). The 
Lipkin-Nogami correction has been applied for ^^"Gd and ^^^Hg. 



serves the particle number. 

The ^^^Hg dynamical moments of inertia arc plotted in 
the upper panel of Fig. 1 within the three mentioned for- 
malisms together with their experimental counterparts. 
As was already the case for -'^"''"Gd, the HF calculations 
fail to reproduce the experimental values as well as the 
behavior of this moment as a function of the angular ve- 
locity since it remains almost constant over the whole 
rotational band. By adding pairing correlations within 
the HFB+LN formalism, it is possible to reproduce the 
experimental values at low spins. However, using a sim- 
ple seniority force as done here, one can only reproduce 
the increasing trend of the experimental moment of in- 
ertia and our theoretical value becomes to high at the 
upper bound of the rotational band. In view of some 
residual instabilities of the solution due to the presence 
of single-particle states in the vicinity of the cutoff en- 
ergy defining the configuration space where pairing cor- 
relations are treated, we have introduced for this nucleus 
a cutoff factor of the Fermi-function type. 

Similarly to what we have done for the ^^"Gd nucleus, 
we have attempted to model the collective effects of pair- 
ing correlations in this nucleus by performing HF+V cal- 
culations. As seen in Fig. 1, the calculated HF+V dy- 
namical moments of inertia nicely match the HFB+LN 
ones, even better than for ^^"Gd. The axial deformation 
of the ^"^Hg nuclear states in the rotational band calcu- 
lated within these two formalisms differ by less than 3 
percent with an axial quadrupole moment mean value of 
Q20 — 44 barns. 

In view of the rather low expectation values of the 
Kelvin circulation in the HF regime (see Fig. 2) its behav- 
ior as a function of the angular momentum in the HFB 
approach reflects mostly the interplay between the global 
rotation and the pairing correlations. One observes in- 
deed a very interesting parabolic pattern. At low spins 
its low value is due to the fact that the intensity of the 
intrinsic currents is obviously related to the intensity of 
the global rotational currents which they are counteract- 
ing. On the other hand these intrinsic currents are, as we 
have assumed and demonstrated here, related with the 



intensity of pairing correlations. This balance explains 

the appearance of a maximum for Ki^i)] at some "inter- 
mediate" spin value (around / ~ 20 fi). From the plotted 
value of (Ki) for both charge states in Fig. 2, one may 
infer that as it was the case for ^'"^"Gd, the pairing gener- 
ated intrinsic currents are almost entirely of an isoscalar 
character. 

The agreement between HFB and HF+V dynamical 
moments of inertia clearly points into the direction of 
identical current patterns for both theoretical descrip- 
tions. This is indeed so, as shown in Fig. 3. Finally, 
one observes in Fig. 4 that the similarity of the kinetic 
moments of inertia in the two theoretical descriptions 
is reflected in the corresponding curves giving the spin- 
dependence of the angular velocity. 

C. The ground-state deformed band of ^^''No 

Recently, the existence of the rotational ground-state 
band of the heavy nucleus ^^^No has been assessed [30] 
by an observation of the 7-transition energies from spins 
J = 4 ?i up to / = 20 fi, in coincidence with the a-decay 
chain of this nucleus. This high spin structure is the first 
one studied in Z > 100 nuclei. 

It has been shown by several authors [14,31,32] with 
diflferent effective forces and pairing parameterizations 
that the properties of this isotope could be reproduced 
within HFB-type formalisms. The calculations presented 
here have been performed in the spirit of Ref. [14] , that 
is fixing the pairing-force parameterization in order to 
reproduce the first transition energies. As a result, it 
happens that higher transition energies are fairly well re- 
produced. This procedure has been found in [14] to be 
more successful within a pure HFB than a HFB+LN for- 
malism. Thus we will discuss here in terms of intrinsic 
vortical currents the results of our HFB calculations. 

The dynamical moments of inertia obtained within the 
three formalisms HF, HFB and HF+V are plotted in the 
corresponding upper panel of Fig. 1 as functions of the 
angular velocity. As was the case for the two previously 
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studied nuclei, it is still impossible here to reproduce the 
experimental data with pure HF calculations. On the 
contrary, the results of our HFB calculations fit rather 
well in the experimental patterns, and the HF+V results 
fully overlap the HFB ones. A slightly less good agree- 
ment between both theoretical estimates of O*-^-* among 
themselves and with their experimental counterpart is 
obtained as exhibited in Fig. 1. 

In Fig. 2, the HF and HFB estimates of {Ki) for the 
^^''No neutron current distributions are seen to be very 
close. This is clearly related as exhibited in the lower 
part of this figure to the very small amount of pairing 
correlations for this charge state at spin less than 6 h 
(whereas for higher spin values, there are no correlated 
solutions). As for the protons now, the shell effect con- 
tribution to {Ki) (i.e., the HF value) is quite significant, 
actually of the same order of magnitude as the pairing 
contribution (i.e., the difference between the HFB and 
the HF value) . As a very specific feature of the HF plus 
vorticity calculations in this nucleus (versus those above 
reported in ^'"'"Gd and ^^^Hg) one notices that a purely 
isoscalar constraint on (i^i) is not effective in reproduc- 
ing the HFB results. Actually, whereas it is tempting 
to attribute such an isovector behavior for a Z ~ 100 
nucleus (as opposed to Z ~ 80 or less nuclei) to the in- 
creasing relative importance of the Coulomb interaction 
(scaled as Z"^), the observed behavior of 5\{Ki)\ (differ- 
ence between HFB and HF-I-V values of K-ftTi)!) which 
is positive for protons and negative for neutrons is yet 
to be understood (and should be assessed by systematic 
calculations). 

In Fig. 3 intrinsic currents corresponding to two very 
different rotational regimes {1 = 2 and 20 h) are dis- 
played for ^^^No. At low spins the neutron currents are 
small yet roughly consistent with a tangential pattern at 
the surface. This is due to the existence of a low pairing 
correlations regime. At high spin of course the HF and 
HFB neutron current patterns are coinciding and hardly 
tangential at the surface. This situation is quite at vari- 
ance with what is obtained for the protons at large spins, 
where the prevalence of large pairing correlations entails 
the existence of rather well marked tangential counter- 
rotating currents. 

Finally, the behavior of the angular velocity of both 
HFB and HF-I-V calculations on the one hand and of the 
HF calculations on the other hand, as a function of the 
total angular momentum is found quite similar here (see 
Fig. 4) to what was found before for ^^°Gd and ^^^Hg. 

III. CONCLUSIONS 

In the above, we have examined in a Routhian ap- 
proach, the equivalence of the dynamical effects yielded 
by pairing correlations a la Bogoliubov with those ob- 
tained when a specific ansatz for the collective flow is im- 
posed onto a purely Hartree-Fock type of microscopic cal- 



culations. This has been achieved in rather different con- 
texts as far as the nuclcon number [A ranging from 150 
to 250) , the deformation (considering normally deformed 
bands as well as superdeformed bands), the range of spin 
values as well as the pairing correlation contents. In all 
cases however, we have found rather convincingly that 
constraining Hartree-Fock Routhian solutions to have the 
currents obtained by HFB Routhian calculations at the 
same spin values, was a sound substitute to these HFB 
calculations themselves. 

For that purpose, the so-called Kelvin circulation op- 
erator has appeared as an efficient measure of the in- 
trinsic currents. Since this operator is well suited to 
the collective dynamics dubbed as S-type ellipsoids by 
Chandrasekhar (implying a collective field which is linear 
in the coordinates), the actual HFB dynamical pattern 
should not deviate very much from this simple model 
pattern. That has been actually assessed by a mere in- 
spection of current lines obtained in our HFB calculations 
in cases wh(;re at large enough spin values strong pairing 
correlations where present. Pairing correlations appear 
indeed as generating tangential counter-rotating currents 
in the intrinsic frame. This equivalence has been actually 
confirmed by a convergent array of evidences: 

i) the similarity of dynamical moments of inertia, 

ii) the congruence of kinetic moment of inertia, or 
equivalently (at the same value of the angular mo- 
mentum) of the associated angular velocity, 

iii) the correspondence between currents obtained in 
both approaches. 

As for the currents, it has appeared that shell effects 
where visible in situations where either the collective in- 
trinsic flow was too small (at low angular velocity for 
the global rotation) or when pairing correlations where 
vanishing. These shell effects generating density inho- 
mogcneitics are responsible for erratic current patterns 
(sometimes roughly similar to localized vortices) which 
are in general neither counter-rotating nor tangential. 

As above disciissed, upon increasing the spin value one 
experiences a competition between two processes. The 
first one is an increase of the intrinsic currents, a reaction 
effect which should be (and actually is, all things kept 
equal) increasing with the global rotation angular veloc- 
ity. This type of currents decreases, as well known, the 
kinetic energy at a given angular velocity. The second one 
is the Mottelson-Valatin rotation-generated reduction of 
the pairing correlations on. The competition arises from 
the fact that this intrinsic flow should be decreasing with 
the decrease of the pairing correlations. One may thus 
easily understand the maximum obtained for the abso- 
lute value of the expectation value of the Kelvin circula- 
tion obtained for the ^®^Hg nucleus at intermediate spin 
values. 

In our calculations an isoscalar constraint on the 
Kelvin circulation operator has been used. Whereas it 
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has been deemed to be adequate for the two lighter nu- 
clei under study, this has boon hardly the case for the 
heavier one. Confirming and, if so, understanding this A- 
dependence or more likely Z-dependence, would require 
a rather systematic study which has not been performed 
here. 

The present study is now extended in two directions. 

At high spins the decrease of pairing correlations has 
often been extrapolated towards the appearance of a 
phase transition (associated thus with a complete dis- 
appearance of these correlations). It is clearly not very 
appropriate to quantitatively assess this feature in tradi- 
tional approaches because one reaches there a low pairing 
regime where the Bogoliubov particle number symmetry- 
breaking approach is not very adequate there as well 
known. We are currently undertaking [33] Routhian cal- 
culations which explicitly conserve the particle number 
according to the general lines proposed in Ref . [28] . 

Even though we are representing here in a rather sim- 
ple way the dynamical effects of pairing correlations at 
finite angular momentum as obtained in a Bogoliubov 
approach, we are still needing such lengthy calculations 
to yield the correct amount of intrinsic currents (through 
the expectation value of the Kelvin circulation operator). 
The challenge is to try to correlate for a given nucleus 
(i.e., for given nucleonic numbers, given deformations and 
given single particle level densities near the chemical po- 
tential at zero spin or equivalently the amount of pairing 
correlations at zero spin) the Kelvin circulation expec- 
tation value as a function of the angular velocity of the 
global rotation. This is currently achieved through a sim- 
ple model approach inspired from the Mottelson-Valatin 
rotation anti-pairing effect and will be presented soon in 
an other paper [34] . 
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